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(-H 1 Abstract. A polynomial of the form x a — p(x), where the degree of p is less 

' than the total degree of x a , is said to be least deviation from zero if it has the 

smallest uniform norm among all such polynomials. We study polynomials of 
least deviation from zero over the unit ball, the unit sphere and the standard 
simplex. For d = 3, extremal polynomial for (x\x%x$) k on the ball and the 
sphere is found for k = 2 and 4. For d > 3, a family of polynomials of the form 
(xi ■ ■ ■ x d ) 2 — p(x) is explicit given and proved to be the least deviation from 
zero for d = 3, 4, 5, and it is conjectured to be the least deviation for all d. 



1. Introduction 



Let 11^ denote the space of polynomials of degree at most n in d variables and 
we write II n = n*. For d = 1, it is well-known that the 2 1 "™ multiple of the 
Chebyshev polynomial of the first kind 

T n (x) = cosn(arccosx) = 2 n ~ 1 x n + q{x), q e II„_i, 

is the monic polynomial of least deviation from zero in II n in the space C[— X, 1]; 
that is, 



O . pen„_ 



inf \\x»-p(x)\\ci-i,i] = 2 1 -"||T n || c[ _ M] = 2 1 "". 



Equivalently, we say that x n — 2 1 ~ n T n is the best approximation to x n in C\— 1, 1]. 
Let fi be a region in M d . For / € C(O), the best approximation of / from 11^ in 
i^h | the uniform norm is the quantity 



(1-1) #»(/;«)= in| ||/-p||c(fi), 



, where |j/||c(o) — niax xe fj We call p* an extremal polynomial for / if 

E„(/;0) = ||/ -p*\\ C (ay For x = (x 1 ,...,x d ) € M d and a = (ai,...,a d ) € Ng, 
we define the monomial d . The degree of the monomial x a is 

\a\ = a± + . . . + a c i- If p*{x) is an extremal polynomial for the monomial x a , 
we call x a —p*(x) the polynomial of least deviation from zero. For f2 being a region 
in M. d , polynomials of least deviation are known only in the case that SI is a cube. 
We are interested in the case of the unit ball B d = {x : \\x\\ < 1}, where is the 
usual Euclidean norm of x, the unit sphere S^ 1 — {x : \\x\\ = 1} and the standard 
simplex T d = {x : X\ > 0, . . . , Xd > 0, 1 — x\ — . . . — x d > 0}. 

For d = 2, the least deviation of x n y m from II^ +m _ 1 in the space C{Tl) has been 
studied for B 2 and T 2 (see, for example, |Gearhart, 73| , |Reimer, 77| , |Newman and Xu,93| , 
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|Bojanov~ Haussm ann and Nikolov, 01| ). For d > 2, the only case known is xi ■ ■ ■ Xd 
on B d and S d -\ which is a polynomial of least deviation by itself. This is shown 
recently in |Andreev and Yudin, 01| , 

(1.2) inf \\x 1 ---x d -p{x)\\ c(B d ) = inf \\x x ■ ■ ■ x d - p{x)\\ C (s<i-i) 

= Wxf-XdWc^d-i) =d- d/2 . 

In other word, the best approximation of x\ ■ ■ ■ Xd from is the zero polynomial. 
Finding polynomials of least deviation on these regions appear to be a difficult 
problem. Only a handful of explicit non-trivial examples of extremal polynomials 
for d > 3 are known in the literature. 

In the present paper, we study the least deviation from zero for monomials of 
lower degrees. We found extremal polynomials for (X1X2X3) 2 and (X1X2X3) on 
B 3 and S 2 and a family of extremal polynomials for x\ ■ ■ ■ x 2 d on B d and S d ~ 1 , 
which are derived from the extremal polynomials for X1X2X3 and (X1X2X3) on T 3 
and xi ■ ■ ■ Xd on T d , respectively. We give an explicit construction of this family 
of polynomials and conjecture that they are the least deviation polynomials. The 
conjecture is proved for d = 3,4,5. The result provides, we believe, the first non- 
trivial example of polynomials of least deviation on these domains. For example, 
we have 

inf \\xlx%xl-p(x)\\ C (B3) = inf \\x\x\x\ - p(x)\\ c{S 3 ) = 72~ 2 

and the minimum is attained by the extremal polynomial R^^x^y) defined by 
R 3 (x 1 ,x 2 ,x 3 ) = 12x\x\x\ - A{x\ + x\ + xf) + A{x\ + x* + x^) 2 + 1. 

The least deviation, 72~ 2 , is surprisingly small in view of the value 3 3 / 2 for X1X2X3. 

Our proof is based on a general result for the Chebyshev approximation in 
|Rivlin and Shapiro, 61| , in which the best approximation element is characterized 
in terms of extremal signature. The most difficult part, however, is to identify a 
correct extremal polynomial. There is no general method for this purpose. We 
relied heavily on the computer algebra system Mathematica to test and verify con- 
jectures. In retrospect, the explicit construction is natural and rather suggestive. 
For example, Rs(x) agrees with the Chebyshev polynomial T2(x) on the three edges 
of the face of T 3 defined by x\ + x 2 + X3 = 1. The result allows first glimpse of 
what an extremal polynomial in more than two variables may look like. 

The paper is organized as follows. In the following section we recall the theoretic 
background needed to prove our result. The results for d = 3 are discussed in 
Section 3 and those for d > 3 are in Section 4. 

2. Extremal signature and best approximation 

We recall the characterization of the extremal polynomials in terms of the ex- 
tremal signature. The study in |Rivlin and Shapiro, 61| is given in the general 
setting of approximation from a finite dimensional subspace of C(Q) on a compact 
Hausdorff space f2. We shall restrict the statement to our setting. 

Let f2 be an infinite compact set in M. d . A signature a on the set fl is a function 
with finite support, whose nonzero values are either +1 or —1. A signature a is 
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called extremal with respect to 11^ if there exists a subset S in the support of a 
and positive numbers A„, v £ S, such that 

X v a(v)p(v) = 0, for all p in 11^. 

v&S 

Let r > be a fixed number. For each p £ , we denote by S r (p; f) the set 

Srfe f) = {x&n:\f(x)-p(x)\ = r}. 

If r = ||/ — pj|c?(£7) , S r (j>', f) is the set of extremal points of / — p and we denote it 
hy S(p;f). 

The characterization of the best approximation of / from 11^ is given by the 
following theorem in |Rivlin and Shapiro, 61| : 

Theorem 2.1. A polynomial p* satisfies \\f— P*\\c(Q) — E n (f; CI) if and only 

if there exists an extremal signature a with support in S(p*]f) such that cr(v) — 
sign(/-p*)(v) for alive S(p*;f). 

The sufficient part of the theorem provides a method to verify if a polynomial p* 
is extremal. One needs, however, to know the extremal polynomial in advance, as 
the extremal signature is supported on the set S(p*; f) which depends on p* . The 
sufficient part of the theorem can be extended to the signature support on S r (jp; f), 
in which r is not necessarily ||/ — p||c(fi)- We will use this slightly extended version, 
which we state in the following. A simple proof is included for completeness; sec 
|Rivlin and Shapiro, 61| for more details. 

Theorem 2.2. Suppose there exists a polynomial p* £ 11^ and an extremal signa- 
ture a supported on S r (j>*; f). Then E n (f; Q) > r. 

Proof. We can normalize the measure A M for the extremal signature so that it is 
a probability measure; that is, YlveS (p*-f)^ v = ^et ^( r ) = S r (p*if) m this 
proof. Since ^vP{v) = for any polynomial p 6 11^, we have 

\\f(x)~p(x)\\ cm > £ WV)-P(V)\ 
vGS(r) 

> | ^ ^vO- v f(v) - ^ X v a v p(v) 
= | Y K<rvf(v)\ = \ Y X vO- v (f(v)-p*(v)) 
= ^\f(v)-p*(v)\=r ]T \ v =r. 

where we have used the fact that f(v) —p*(v) = o~ v r for v S S(r). □ 

The extension allows us to apply the result to the situation where a good can- 
didate for p* is identified but the norm of / — p* is hard to determine. This is 
precisely our case in Section 4. 

Our construction is motivated by the recent study in |Andreev and Yudin, 01| , in 
which it is shown that if / is invariant under a finite group G (that is, f{xg) = f{x) 
for all g £ G), then the best approximation E n (f]S d ~ l ) is attained at G invari- 
ant polynomials. (This result appeared early in |Ganzburg and Pichugov, 81] as 
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pointed out by a referee.) More precise, we state the result in |Andreev and Yudin, 01| 
as follows: 

Proposition 2.3. Let G be a subgroup of the rotation group 0(d) and let GH d 
denote the polynomials in IT^ that are invariant under G. If f is invariant under 
G, then 

inf \\f(x)-p(x)\\ c , S d-i ) = inf \\f(x)-p(x)\\ c(sd -i ) . 

Using this fact, the best approximation of several invariant functions are given 
in |Andreev and Yudin, 0l| , including the case x\ - ■ -Xd in <|1. 2JI (invariant under 
the symmetric group). The proof in |Andreev and Yudin, 01 can be applied to any 
region Q and / that are invariant under a finite group G. In particular, if / is 
invariant under a subgroup G of the symmetric group Sd(T d ) of the simplex T d , 
then an extremal polynomial of / can be taken as a G-invariant polynomial. 

If / is even in each of its variables, then / is invariant under the sign changes 
of each variable (invariant under the group Zj); the extremal polynomial can be 
taken as a polynomial even in each of its variables. Furthermore, instead of B d or 
S x , we can work with T d and T d_1 in this case. In fact, the following general 
proposition holds: 

Proposition 2.4. Let a S Nq and write 2a — (2a%, . . . , 2ctd) and \a\ = n. 
If p*(x) is an extremal polynomial for E n (x a ;T d ) then p*(x\ , . . . , x^) is an ex- 
tremal polynomial for E n (x 2a ; B d ) ; conversely, if q* is an extremal polynomial for 
E n (x 2a ; B d ) in the form q*(x) — p*(x 2 , . . . , x 2 ,), then p*(x) is an extremal polyno- 
mial for E n (x a ; T d ). Furthermore, let f a (x\, . . . , Xd-i) = x" 1 ■ ■ ■ a^l^ 1 (1 — Xi — 
. . . — Xd-i) ad ; then the above conclusion holds for E n (f a ; T d ~ 1 ) and E n (x 2a ; S' d_1 ). 

We note that f a {x\, ■ ■ ■ i^-i) = x 2a on S' d ~ 1 . The proposition follows easily 
from the fact that x i— > (x\ , . . . , x 2 d ) is one-to-one from T d to = {x E B d : xi > 
0}, and the map also induces a one-to-one mapping from 11^ to GH d n with G = Z d > 
that is, the subspace of polynomials that are even in each of its variables. For d = 2 
the proposition has been used in |Bojanov, Haussmann and Nikolov, 01| . The cor- 
respondence between polynomials on these domains also works for other problems 
involving polynomials, such as orthogonal polynomials and cubature formulae, see 
for example |Xu, 98| . 

3. Least deviation from zero for d = 3 

We consider best approximation to the monomials X1X2X3 and (X1X2X3) 2 in this 
section. The main task is to identify an extremal polynomial. The results in the 
previous section provide some guidance, but there is no general method for this 
purpose. Our first example, Rs(x) given below, was found after many attempts. 
See the comments after the proof. 

Theorem 3.1. Define the polynomial Rs(x) by 

i? 3 (a;) = 72xix 2 x 3 — 4(x\ +x 2 +x 3 ) +4(xi +x 2 +x 3 ) 2 — 8(xix 2 + x 2 x 3 + X1X3) + 1. 
Then 72 _1 i?3(a;) is a polynomial of least deviation from zero and 

E 2 { Xl x 2 x 3 ;T 3 ) = £2(2:12:2(1 -xx- x 2 )-T 2 ) = 72- 1 ||i? 3 || c(T 3 ) = 72" 1 
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Furthermore, 72 R 3 (x\, x 2 , x 2 ) is a polynomial of least deviation from zero and 
E 5 (x 2 x 2 x 2 ;B 3 ) = E 5 (x 2 x 2 x 2 -S 2 ) = 72~ 2 \\R 3 {x\,xlx 2 )\\ c{B3) = 72~ 2 . 

Proof. By Proposition ^. 4l we only need to work with the simplex. It is easy to verify 
that ,R 3 (0,0,0) = 1 and R 3 {\/2, 1/2, 0) = -1. Solving the equations diR 3 (x) = 
0, i = 1,2,3, shows that R 3 has 4 critical points inside T 3 but none of them is 
maximum or minimum, since the values of |i?3(x)| at these points are less than 1. 
Thus, I J?3 (ac) I attains its maximum on the boundary of T 3 . It is easy to verify that 
the polynomial R 3 (x) satisfies 

R 3 (x, y, 0) = R 3 {x, 0, y) = R 3 (0, x, y) = (1 - 2x) 2 + (1 - 2y) 2 - 1 

which is bounded by 1 in absolute value. Hence, we only need to show that |i?3(x)| 
is bounded by one on the face of T 3 defined by X \ + x 2 + X3 — 1; that is, we need to 
show that U 3 (x\, x<z) = R 3 (x\, x 2l 1 — x\ — x%) is bounded by 1 in absolute value on 
T 2 . Taking derivatives of U 3 (xi, x 2 ) and solving for the critical points shows that 
it has 4 critical points inside T 2 , of which only the point (1/3, 1/3) is a maximal, 
E/(l/3,l/3) = 1. Furthermore, it is easy to verify that 

U 3 (x,0) = U 3 (0, x) = U 3 (x,l-x) = T 2 (x); 

that is, it agrees with Chebyshev polynomial of degree 2 on the boundary of the 
triangle. Hence, |[/ 3 (a;)| < 1. Furthermore, the above analysis also shows that 

S+ := {x : \R 3 (x) = 1} = {(0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), (1/3, 1/3, 1/3)} 
5_ :={x:\R 3 (x)=-l}= {(1/2, 1/2,0), (1/2, 0,1/2), (0,1/2, 1/2)}. 

Let <j(v) = 1 on S+ \ {(0,0,0)} and a{v) — —1 on 5_. We show that a is an 
extremal signature. Define L\f and L 2 f by 

Lif = ^ Q, ~ \) + 0, 0) + f(0, 1, 0) + /(0, 0, 1)) 



and 



W = ^(/('i,i,o) + /(i,o,i)+/(o,i > i 



Then Lf := L\f — L 2 f satisfies Lf = 0, / E H 2 . Thus, u is an extremal signature 
for x\x 2 x 3 in C(T 3 ). Furthermore, the support sets S+ and 5_ of a are on the face 
of T 3 , which is identified with T 2 . This shows that a is also an extremal signature 
for Xl x 2 (\ - xi -x 2 ) in C(T 2 ). □ 

Let us mention a connection between cubature formulae and the extremal signa- 
ture for i?3(x). A cubature formula is a linear combination of function evaluations 
that gives an approximation to an integral ( |Stroud, 71| ). Let d/j, be a positive 
measure on 17 C M. d . If 

JV 



f f(xW = ^\ k f( Xk ), feut 



k=l 

rd 



and there is at least one / € n„ +1 such that the equality fails, then the cubature 
formula is said to be of degree n. It is called positive, if all are positive numbers. 
For the extremal signature for R 3 (x), it is easy to verify that both L\f and L 2 f are 
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cubature formulae of degree 2 for dx on the set S 2 = {x e T 3 : x\ + X2 + x 3 = 1}; 
that is, 

/ f(x)dx = L 1 f = L 2 f, for all/ in n|. 

Since we identify S 2 with T 2 , one can write L\ and L 2 as linear combinations of 
function evaluations for functions of two variables. Thus, the extremal signature is 
given by the difference of two positive cubature formulae. 

During our search for R 3 (x), we found U 3 (x) first. In retrospect, the formula of 
U 3 (x), which can be written as 

U 3 (x\,x 2 ) = 72xia;2(l - xi - x 2 ) - 3 + 4(x 2 + x\ + (1 - xi - x 2 ) 2 ), 

is quite natural since it agrees with the Chebyshev polynomials of degree 2 on the 
boundary of T 2 . This also suggests the possibility that other monomials may also 
have extremal polynomials that agree with Chebyshev polynomials on the bound- 
ary of T 3 . For example, for (xix 2 x 3 ) n , one may look for a polynomial that agrees 
with Chebyshev polynomials of n-th degree on the boundary of the simplex. One 
example of such a polynomial is T n (R 3 (xj), where T n (t) denotes the Chebyshev 
polynomial of degree n. Although this function is not a polynomial of least devia- 
tion, it helps us to find a solution for the monomial (xix 2 x 3 ) 2 . 

Theorem 3.2. Define polynomial by 

R 5 (x) =27 2 b (x lX2 x 3 ) 2 - 1 + 2(2;! + x 2 + x 3 ) - 2( Xl + x 2 + x 3 ) 2 

+ 2 [l - 4(xi + x 2 + x 3 ) + A{x\ +xj+x 2 3 )] 2 

- 27xix 2 x 3 [(32/9 - 2a + b)(xi + x 2 + x 3 ) 2 + 6a(xix 2 + xix 3 + x 2 x 3 )] . 
Then 21~ 2 b~ 1 R${x) is a polynomial of least deviation from zero, 

E 5 (x 2 x 2 x 2 ;T 3 ) = E 5 (x 2 x 2 (l - Xl - x 2 ) 2 ;T 2 ) = 27- 2 b- 1 \\R 5 \\ c{T3) = 27~ 2 b-\ 
where the constant a, b and the reciprocal of the least deviation is given by 

a = 28.5926243, 6 = 21.8935834, 27 2 6 = 15960.4223. 
Furthermore, 72~ 4 b~ 2 R 3 {x\ polynomial of least deviation from zero, 

E 5 (x 2 x 2 x 2 ;B 3 ) = E 5 (x 2 x 2 x 2 ;S 2 ) - 72- 4 &- 2 ||i? 5 (x 2 , x 2 , x 2 ) \\ C(B3) - 27- 4 &- 2 . 

Just like the case of R 3 (x), the proof amounts to showing that |-R5(x)| < 1 on 
T 3 and there exists an extremal signature. It is not difficult once the formula of R5 
is identified. We first give an account on how R 5 is discovered. 

Following the construction of R 3 (x), we look for a polynomial in the form of 

U 6 (x, y) =21xy{l - x - y) [27bxy{\ - x - y) + 3a(a; 2 + y 2 + (1 - x - y) 2 ) - c] 

+ 2 [-3 + 4(x 2 + y 2 + (1 - x - y) 2 )} 2 - 1 

that will be a polynomial of least deviation on T 2 with leading monomial x 2 y 2 (l — 
x — y) 2 . Note that the Chebyshev polynomial of degree 2 is T 2 (t) = 2t 2 — 1 and 
T 2 (2t-1) = -3 + 4(< 2 + (l-t) 2 )). The form of U 5 is chosen so that on the boundary 
of T 2 it satisfies 

U 5 (x, 0) = U 5 (0, x) = U 5 (x, l-x) = T 2 (T 2 (2x - 1)) = T 4 (2x - 1). 
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We then choose c = 2/9 + a + 6so that U 5 (l/3, 1/3) = 1. It follows that 1-U 5 {x,x) 
can be factored as 

1 - U 5 (x,x) = x(l - 2x)(l - 3x) 2 (64 - 54ax + 27bx + 162bx 2 ). 

We need to choose a and b so that the last factor is positive for < x < 1/2. One 
choice is to make this factor 2b(9x + d) 2 . This leads to a = 16(3 — Ad) /(3d 2 ) and 
b = 32/d 2 . At this point, it becomes apparent that there need to be more points on 
which U$(x, y) = — 1 inside T 2 . We therefore solve the equations Uc,{x, x) = — 1 and 
C/5 (x, x) = 0. This leads to d = —1.208972894, which gives the values for a and b in 
the theorem. It turns out that this choice does work out and ^(x, y)| < 1 on T 2 . 
The final step is to identify the formula of i?s(xi, X2, x$) from that of U§(xi,X2) 
with the requirement that R^(xi,X2, 1 — X\ — x^) = U$(x\,X2) and \R^(x,y)\ < 1 
on T 3 . This step is not trivial since an additional multiple of (.ti + a; 2 + xs) k to 
any term in R$ (x) does not change the value of the polynomial on the face of T 3 
defined by 23 — 1 — x\ — X2- The polynomial {/^(xi, X2) is an extremal polynomial 
on the triangle T 2 that agrees with the Chebyshev polynomials of degree 4 on the 
three boundary lines of T 2 , its graph is depicted below. 




The polynomial C/3 

Let us point out that there does not seem to be a closed form for the values of 
a and b. In fact, the value of d in the above paragraph is one of the real roots of 
the following polynomial: 

- 612220032 - 1365527808 1 - 835528041 1 2 - 101556504 1 3 



23270976 1 + 26037504 1 5 + 7670016 1 



929280 i 7 + 41984 1 8 . 



This polynomial has 4 real roots and 4 complex roots, and it cannot be factored 
over the integers. 

We now give a formal proof of Theorem 13. 21 
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Proof of Theorem l?Ol First of all, we need to show that |i?5(a;)| < 1 for x G T 3 . 
Solving diRz(x) = 0, i = 1,2,3, numerically for critical points shows that |i£s(a;)| 
attains its maximum on the boundary of T 3 . Furthermore, 

R5(x,y,0) = R 5 (x,0,y) = R 5 (0,x,y) 

= -1 + 2{x + y) - 2{x + y) 2 + 2(1 - 4(x + y) + 4(x 2 + y 2 )) 2 , 

and the polynomial has no critical point inside T 2 . Consequently, the maximum of 
|i2s(a;)| is attained on the face of T 3 defined by x\ + x 2 + x 3 = 1. In other words, 
we only need to show that \U~5(xi, x 2 )\ < 1 on T 2 . Again this can be proved by 
solving diU§(xi, x 2 ) = 0, i = 1,2, and the maximum is attained on the boundary. 
This proves that |i?5(x)| < 1 on T 3 and it also gives the set S+ = {x : |i?5(x)| = 1} 
and the set 6>_ = {x : |-Rs(£)| = — !}■ Let S3 be the symmetric group of three 
elements. For a — (01,02,03) € M 3 we define or := (a Tl , a T2 , a T3 ) , r e S3 and 
(o)g := {or : r 6 S3}. Then 

S+ = {(1/3, 1/3, 1/3), (0, 0, 0), (1, 0, 0) G , (1/2, 1/2, 0) G , (t u h, 1 - 2h) G } 
S_ = {((2 - v/2)/4, (2 + V2)/4, 0) G , (t 2 , t 2 , 1 - 2t 2 ) G }, 

where h = 0.4588164122 and t 2 = 0.1343303216. We consider the signature a 
defined by cr{v) = 1, v £ S + \ {(0, 0, 0)} and a(v) = —1, v G S_. To show that a is 
an extremal signature, we define Lf by 

Lf =co/(l/3, 1/3, 1/3) + ci /((1, 0, 0)r) 

r 

+ c 2 X] /((1/2, 1/2, 0)r) + c 3 X] f((ti,ti, 1 - 2ti)r) 

r r 

- c 4 X] /(((2 - V2)/4, (2 + V2)/4, 0)r) - c 5 ^ f((t 2 , t 2 , 1 - 2t 2 )r), 

T T 

where the sum is taken over all distinct permutations of the base point and the 
coefficients are given by 

c = 0.0997251873, ci = 0.0097228135, c 2 = 0.0621246411, 
c 3 = 0.0243979796, c 4 = 0.0615774830, c 5 = 0.1178707075. 

Then Lf = for all / G n 3 , which shows that cr is an extremal signature. This 
completes the proof of Theorem 13.11 □ 

The linear functional Lf given above is evidently a sum of two linear functionals 
with positive coefficients. Unlike the case of R3, however, the two linear functionals 
are not cubature formulas of degree 5 with respect to the Lebesgue measure. 

The two cases solved in this section appear to indicate a surprisingly complicated 
picture for the best approximation of monomials in three variables, and the picture 
is remarkably different from that of one and two variables. We make two remarks 
in this regard. 

Remark 3.1. One surprising fact of Theorem 13.21 is that the least deviation is not 
given by a reciprocal of an integer. This indicates a major difference between the 
case of three variables and that of one and two variables. In the case of one variable, 
the polynomial of least deviation from zero is the classical Chebyshev polynomial, 
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for which the least deviation of x n to II„_i in C[— 1, 1] is 2 1 n . In the case of two 
variables, we know for example 

E n (x k y n - k ;B 2 )= inf \\x k y n - k - P (x)\\ C (b*) = 2 1 "". 

For three variables, however, we do not know if the least deviation of x a ton,, , 
could be represented by a simple formula that depends only on the total degree of 
the monomial. The result in this section seems to indicate that such a formula does 
not exist. 

Remark 3.2. The values of the least deviation in Theorems 13.11 and 13 . 21 are surpris- 
ingly small. Let us examine the case of the unit ball. We know 

(3.1) E n (x k x%- k ;B 3 ) = E n (x k x^ k ;B 3 ) = E n {x k x n ^ k - B 3 ) = 2^ n . 

This follows from the fact that an extremal polynomial p* for x™x 2 must be even 
in 2:3 since x™x 2 is invariant under the group 1 2 applied on the third variable. Let 
p* be so chosen; then 

\K- m x^-p*(x 1 ,x 2 ,x 3 )\\ c(B3) > \\x?- m x? - p*(x 1 ,x 2 ,y/l- x\ - xl)\\ c[B2) 

> inf \\x?r m x? - P (x u x 2 )\\ c{ b^ = 

Furthermore, the equality holds since an extremal polynomial for Xi~ m x 2 n on B 2 
can also serve as an extremal polynomial on B 3 . Below is a list of other cases that 
we know on the unit ball: 

E n (x lX2 x 3] B 3 ) = 3~ 3 / 2 , E n {x\x\x\\B 3 ) = 2~ 6 • T 2 , 
E n {x\x\x\; B 3 ) = 0.5340799374- 2" 12 • 3" 12 , 

where we rewrite the value of the third one, which is given in Theorem 13.21 for 
easier comparision. The value of E n (xlx 2 xl; B 3 ) appears to be strikingly small. 
For other degree 6 monomials given in (33) , the value of the best approximation 
is only 2~ 5 . Also note the fast decrease shown in these three values. 

4. Least deviation from zero for d > 3 

We consider the best approximation to (x± ■ ■ ■ Xd) 2 on B d or S' <i_1 , and the best 
approximation to Xf-Xd on T d or T' 1 ^ 1 in this section. The extremal polyno- 
mial can be taken as symmetric polynomials by Proposition 12. 31 It is well known 
that every symmetric polynomial can be written in terms of elementary symmetric 
polynomials ( jMacdonald, 95| ). 

The elementary symmetric polynomials of degree k in variables x\,x 2 ,... ,xn 
are defined by 

e k (x) = 2J x tl x i2 ■ ■ ■ x ik , l<k<N. 

l<i 1 <...<i k <N 

In particular, ei(xi, . . . , xn) = x\ H + xn and ed{x\, . . . , xn) = x\ - ■ ■ xn- As 

it is often the case with the symmetric functions, we assume that N is sufficiently 
large and do not write the dependence of e& on the number of variables. We will 
use the notation l k = (1, 1, . . . , 1) e M fc . 
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Definition 4.1. Using elementary symmetric functions, define T 3 (x) by 

T 3 {x) = 72e 3 (x) - 4ei(ar) + 4e 2 (x) - 8e 2 (x) + 1 

and Tk(x) for k > 3 by the recursive formula 

T k (x) = r k e k (x) - T fe _i(a:), 

where the constant r k is determined by r k = k k [Tk-i(k~ 1 l k ) + 1]. 

Note that fc _1 l fe = (fc -1 , . . . , k^ 1 ) e R fc ; we use the evaluation of T k -i, as a 
function of R k , at this point in the definition of r k . Clearly r k is uniquely deter- 
mined. For x G R d , the function Td(x) will serve as extremal polynomials. In 
particular, the polynomial T 3 (x) for x e M 3 is the same as R 3 (x) in the previous 
section. For x <G M 4 , the explicit formula of T 4 is given by 

Ti(x) = 896x1x2X3x4 — 72(xiX2X3 + X1X2X4 + X1X3X4 + X2X3X4) 
+ 4(xi + x 2 + x 3 + x 4 ) - 4(xi + x 2 + x 3 + x 4 ) 2 
+ 8(xi.x 2 + X1X3 + X1X4 + x 2 x 3 + x 2 x 4 + X3X4) - 1. 

The value of rd is of particular importance. It can be computed using the following 
formula. 

Lemma 4.2. For d > 3, 

r d = dVfc" ( f) [(-l) fe (9fc 2 - 32k + 24) + fc 2 ] . 

In particular, r 3 = 72, r 4 = 896, r 5 = 14400, and r 6 = 283392. 

We defer the proof to the end of the section and continue to state our main result 
of this section. 

Theorem 4.3. For d > 3, on the d-dimensional simplex 

£ d _i(xi ■■■x d ; T d ) = £ d _i(xi • • • x d _i(l - x x x A - X );T d - x ) > 

and the equality holds for d — 3,4,5 with r^ 1 Td(x) as a polynomial of least deviation 
from zero. Furthermore, on B d and S 4-1 , 

E 2 d-i(x\ • • • x 2 ; B d ) = E 2d - X {x\ • • • x 2 ; S^ 1 ) > r^ 2 

and the equality holds with r d ~ 2 Td{x\, . . . , x 2 d ) as a polynomial of least deviation from 
zero. 

We believe that the equality still holds for d > 6. In fact, all that is missing is 
to prove that |Td(x)| < 1 for x e T d . We state it as a conjecture. 

Conjecture 4.4. For d> 6, the inequality \\Td(x)\\c( T d<) < 1 holds. In particular, 
the equality in the above two theorems holds for d > 6. 
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Let us point out that there does not appear to exist a closed formula for r d . 
Below is a list of the first few values of r d and their prime factorization: 

r 3 = 72 = 2 3 • 3 2 , r A = 896 = 2 7 • 7, r 5 = 14400 = 2 6 • 3 2 • 5 2 , 

r 6 = 283392 = 2 s • 3 3 • 41, r 7 = 6598144 = 2 9 • 7 2 • 263, 

r 8 = 177373184 = 2 15 • 5413, r 9 = 5406289920 = 2 12 • 5 • 3 4 • 3259, 

no = 184223744000 = 2 14 • 5 3 • 23 • 3911, 

m = 6939874934784 = 2 14 • 3 3 ■ ll 2 ■ 137 • 409. 

A closed formula will have to catch the pattern of the prime numbers presented in 
these formulae, which seems unlikely. We also note that the values of r d appear 
to indicate that the best approximation to monomials becomes increasingly more 
complicated as d increases. See also the Remark 3.1. 

The proof of Theorem 14.31 is split into several propositions. As before, we only 
need to prove the case of polynomials on the simplex. We start with the point set 
at which T d (x) = ±1. 

We will work with sets of points that are invariant under the symmetric group 
S d . For a £ M. d we use the notation (a)c to denote the set of points that consist of 
all distinct permutations of x; that is 

(ai, . . . ,a d ) G = {(an, ■ • - ,a-rj : r = (n, ...,r d )£ S d } 
and we sometimes write ar = (a T1 , . . . , a Td ) for r £ S d . 

Proposition 4.5. Let S+(T d ) andS-(T d ) be the subsets ofT d on which T d (x) = 1 
and T d (x) = —1, respectively. For odd d, 

5+ = |(i,...,i) G , (333,... , 3=2' '°) G >- ■■ > (h b 3,0, .. .,0) G ,(1,0, . ..,0) G | 
is a subset of S+(T d ) and 

S - = {(d=T' ■• •> 3=7 '°) G ' (d=3' • •• > 3=3'°'°'°) G " " > (l> 5'°' • • • >°)g} 
is a subset of S~(T d ). For even d, 

= { (d' " ' " ' i)g ' (rf=2 ' ■ ■ ■ ' 3=2 ' 0> ^) G ' ' ' ' > (2' 2' ^' • • • >^)g} 
is a subset of S+(T d ) and 

^- = {(3 1 T--- I 3=T'0) G ,( I i 3,...,^,0,0,0) G ,...,(l,0,...,0)G} 

is a subset of S~(T d ). Furthermore, all these points are on the face of T d defined 
by the equation X\ + . . . + x d = 1. 

Proof. In the definition of T d , the value of r d is chosen so that T d ( -jtj > • • • j ^qrr ) = 1 • 
All other points in the given set contain at least one zero component. This allows 
us to use induction. For T 3 (x), it is easy to verify that T 3 (x) = 1 if x = (1,0, 0)g 
and x = (1/3,1/3,1/3), and T 3 (x) = -1 if x = (1/2, 1/2, 0) G . The induction is 
based on the formula 

T d (xi, . . .,Xd-i,Q) = -T d _i(xi, . . .,x d -i) 

and similar formulae obtained by a permutation of (xi, . . . , x d —i), which follow 
from the definition of T d and the fact that e d (xi, . . . , x d —i, 0) = 0. □ 
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Proposition 4.6. The signature a, defined by o~(v) = 1 ifvE S+ and cr(v) = — 1 
if v G is an extremal signature of Td. More precisely, define 

+ . . . + f-lf-V-' V> (Q, 1,0, ■ -- ,o) t) + YJ 9 ((1,0,. . .,0)r); 

then hp = /or o/Z p G H c j _ 1 . 

Proof. Since the points in 5+ and iS_ are symmetric with respect to Sd, a moments 
reflection shows that we only need to verify Lg — for symmetric polynomials. 
One basis of symmetric polynomials in II^ j consists of mk(x) = x\ + ■ • • + x\ for 
k = 0, 1, . . . , d- 1. We show Lmj = 0. Let ^ = (1/j, . . . , 0, . . . , 0) G T rf , which 
contains exactly j nonzero entries. Then the sum ~^2 TeSd g(ajr) contains (^) terms 
and mfc(aj) — j(l/j) k - Consequently, for k > 1, 

ta, = d d - fe - (f) (d - l) rf - fc + (d - 2) d - fc + . . . + ( d _\ 

Furthermore, it is easy to see that LI gives the same formula as Lm\ (recall that 
points in S+ and <S_ satisfy X\ + . . . + Xd = 1). Hence, we need to show, changing 
n — k to k, 




k = l,2,...,d-l. 



This is well known and can be proved by induction on d. □ 

By Theorem 3.2, the above proposition has proved that 
inf \\x l ...x d -p{x)\\ c{T d ) > rj 1 . 

In order to show that the equality holds, we only need to prove that |Td(a;)| < 1 for 
x G T d . However, we are able to establish this inequality only for d = 3, 4, 5. 

Proposition 4.7. For d — 3,4, 5, 

\Td{x)\<\\T d \\c ( T«-r ) = l, x G T d . 

Proof. There does not seem to be an easy way of proving this. We use the standard 
method of finding critical points upon solving (diT^(x) / dxi) = 0, 1 < i < d. The 
case d = 3 is in the previous section and the equations can be solved algebraically. 
The cases d = 4 and d = 5 are solved numerically. The details are omitted. Once 
the critical points are found, we can then verify that the inequality |Td(x)| < 1 
holds on these points, which shows that the maximum of |Td(a;)| is attained on the 
boundary of T d . Since Td(xi, . . . , Xd-i, 0) = —Td-i(x\, . . . , Xd-i), by induction, we 
only need to prove that \Td{x\, . . . , Xd—i, 1 — X\ — . . . — x d -i)\ < 1 for x G T d_1 . 
Again, this is done by computing the critical points and evaluating. □ 

Putting the above propositions together, we have completed the proof of Theo- 
rem EH 
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We still need to prove Lemma PT21 First we note that the definition of T d implies 

d 

T d (x) = 2(-l)^ fc r fc e fc (a;) + (-l) d - 3 T 3 (x). 

fc=4 

Proof of Lemma \4-i?\ Setting x = a d+ i = (d+ and using the fact that 

efc(a<j+i) = l) _fe leads to the relation 

(d + l)- {d+1 ^ d+1 = ^(-l)""* + k l ) (d + l)- k r k + (-l) d - 3 T 3 (a d+1 ) + 1. 
Replacing d + 1 by d we can write the above equation as 

J2(- r ) d ~ k (t) d ~ krk = (-^) dT s(d- 1 l d ) + 1 := A d . 

k=i ^ ' 

We want to reverse this relation so that r d = Ylj=4,bjAj- A standard argument 
shows that, with bj given below, this will follow from the combinatoric relation 

j—k 

where Sk,d = 1 if k = d and otherwise. For k = d, J dd = 1 holds trivially. For 
k < d, a change of summation index gives 




d-l-k 



By H4.1fl . the last formula shows that Jk,d — 0. Finally, the definition of T 3 shows 
that 

T 3 (^ 1 l d ) = 72 (f\ d- 3 - 8 (f\ d- 2 + l = d- 2 (9d 2 - 32d + 24), 

which gives the explicit value of A d . Putting these relations together completes the 
proof. □ 

We conjecture that |T<j(a;)| < 1 for all d and the points in S+ and S- are the 
only ones on the face of T d defined by x\ + . . . + x d = 1 on which |2d(a;)| = 1. The 
following fact is helpful for the case d = 5 and could be useful for d > 5. Let 



D d (x) = dot 



1 


1 


1 


X\ 


x 2 


Xd 


x d-3 


x d-3 








x d 


pxT d {x) 


diT 2 (x) . 


■ d d T d (x)_ 



where di — d/dxi. Then D d (x) can be factored completely. We have, for example, 
D 5 (a;) = -64 J| ( Xi - Xj)(-U + 225x 5 ). 

l<i<j"<4 

At the critical points of T 5 , D 5 (x) — and D 5 (xt) — for r £ S5. One could use 
it to confirm the conjecture for d > 5. We did not try hard to push for larger d, 
since the method does not seem to lead to a proof for all d. 
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Our conjecture implies that Td(x) attains its maximum on the boundary of T d . 
Part of this can be proved as follows: Let A denote the Laplacian operator A = df + 
. . Then it is easy to verify that ATd(x) — (— l) d_1 8. In particular, this shows 

that (— 1) 1 Td(a;) is a subharmonic function. Hence, by the maximum principle 
( | John, 82] ) for the subharmonic functions, we can conclude that (— l) d_1 T c ((a;) < 
max x g dT d ( - 1 ) d ~ 1 T d ( x ) . 

Acknowlegement. The author thanks the referees for their careful review. 
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